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Abstract 

We study a generalization of the porous medium equation involving nonlocal terms. More precisely, explicit self- 
similar solutions with compact support generalizing the Barenblatt solutions are constructed. We also present a 
formal argument to get the LP decay of weak solutions of the corresponding Cauchy problem. 

Resume 

Solutions auto-similaires pour une equation des milieux poreux non locale 

Cette note est consacree a l'etude d'une generalisation non locale de l'equation des milieux poreux. Plus precisement, 
on obtient des formules explicites de solutions auto-similaires a support compact qui ressemblent fortement aux solu- 
tions de type Barenblatt. On donne aussi un argument formel qui permet d'obtenir des estimations L p des solutions 
faibles du probleme de Cauchy. 



Version franchise abregee 

Nous considcrons le probleme de Cauchy pour l'equation non locale suivantc 

d t u-V-(uV a - 1 (\u\ m - 1 ))=0, (1) 

avec m > 1, x £ R d , a £ (0, 2), t > 0, a laquelle on ajoute une condition initiale w(0, x) = uo(x). Ici V" est 
un operateur integral singulier generalisant le gradient usuel (/? = 1) et lie au laplacien fractionnaire. 

Le resultat principal de cette note sont des formules explicites de solutions auto-similaires qui se propagent 
a une vitesse fmic. 

Theoreme 1 (Solutions auto-similaires) 

La fonction 



u(t,x) = Ct~ ■*<">-!>+« (^[r 2 - \x\H~ Mm.-i)+o,y^ 



(ainsi que ses translations en x) est une solution auto-similaire de l'equation (1) pour R > quelconque et 
une constante C > convenable. 
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Ensuitc, nous prcscntons un calcul formel qui permet d'obtenir des estimations en norme L p (M. d ) des 
solutions faibles du problemc de Cauchy, en particulier ccllcs construites par Caffarclli ct Vazquez [3] dans 
le cas 771 — 2 et telle que |uo(x)| < Ce~ c ^ pour deux constantes C et c. 

Calcul formel 1 (Asymptotique du probleme de Cauchy) 

Etant donnee une fonction < uq G D L°° (M. d ) , les normes L p (R d ), 1 < p < oo, des solutions faibles 

u tendent vers quand t — > oo avec le taux algebrique suivant 

a + d(m-l)/p 

||«(*)||p < C*(d,a,m,p)||ii j| 1 Q+d(m - 1) for all t > 0. 

1. Introduction 

We study a nonlocal generalization of the porous medium equation 

a t «-v-(«v a - 1 (i«r- 1 )) = o, (i) 

where m > 1, a G (0, 2), x G R d , < > 0, supplemented with an initial condition 

u(0,x) = uo(x). (2) 

The pseudodifferential (vector- valued) operator V 3 in (1) is defined via the Fourier transform as V^u = 
-F~ 1 (i£\£\P~ 1 J-u). This definition is consistent with the usual gradient: V 1 = V; the components of V° 
are the Riesz transforms; moreover we have V • V"" 1 = Vs • V' = — (— A)^, where (— A)"2 denotes the 
fractional Laplace operator: (— A)%u = J 7 ^ 1 (\£\ a Tu) . It can also be defined by real analysis tools as follows 
V Q_1 u = Vl2- a u, where la for /3 G (0,d) is an integral smoothing operator, called the Riesz potential (see 
[12, Ch. V]) 

with some Cp > 0. Note that then V a ~ 1 u(x) = Vh- a {u){x) = (d + a- 2)C2- a J (u(x + z) - u(x)) , z ^+ a dz, 
a G (0,2). 

Eq. (1) can be interpreted as a transport equation of the type dtu = V ■ (itv) for some velocity vectorfield 
v which is a potential; more precisely, v = Vp where p = l2- a (\u\ m )■ It can be interpreted as a nonlocal 
pressure in the case of nonnegative initial data. Then, of course (— A) - 2— p = u| m_1 , see [3], [4] for that 
notation. 

Notice that for a = 2 and nonnegative initial data we recover the Boussinesq equation (m = 2), and the 
usual porous media equation (m > 1): dtu = V • (uV (w™" 1 )), t > 0, x G M d . 

Recently, L. Caffarclli and J. L. Vazquez ([3], [4]) studied equation (1) in the case m = 2. They proved 
the existence of weak solutions for nonnegative bounded integrable initial data with exponential decay at 
infinity. They also treat the case of bounded and compactly supported initial data, which propagate with 
finite speed. It is shown in [3] that self-similar solutions can be constructed by considering an obstacle 
problem for the fractional Laplace operator. In this note, we contribute to those results constructing explicit 
compactly supported self-similar solutions. Moreover, we show a kind of hypercontractivity estimates, i.e. 
the optimal decay in LP of general solutions with uq G i 1 (R d ). 

Results on equation (1) in this note are multidimensional generalizations of those obtained in [1] for 
a model of the dynamics of dislocations in crystals (for the integral of u when d = 1 and m = 2). The 
structure of (1) suggests that it should enjoy the conservation of mass and some comparison properties as 
was shown in [1] . For an analysis of a related nonlocal equation, see [7] . 

Complete proofs of all results announced in this note will be published in [2] . 

2. Self-similar solutions 

The equation (1) has the following scaling property: 
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if u(t,x) is a solution, so is £ u(£t,£ x) 

for each t > and A = Wnv-[j+a ' ^ e f° r nonnegative solutions that arc invariant under that scaling, 
i.e. of the following form 

«M = ^*-,m(J), where A = d(m _\ } + Q , (3) 

for a function $ Q . m : M d — > R + satisfying the following nonlinear and nonlocal equation in R d 

-AV • (y$ a , m ) = V • V - 1 *™- 1 ) . (4) 

Before stating our main result, we recall the definition of weak solutions for (1) introduced in [3] in the case 
m = 2. 

Definition 2.1 (Weak solutions) A function u : (0,T) weak solution of (1) in Qt = 

(0,T)xR d submitted to the initial condition u(0 , x) = u (x) i/u G L^Qt), J 2 -a(M m_1 ) G ^(O, T; W 1 x ' c 1 (R <i )), 
uVJ 2 -a(|u| m_1 ) G L^Qt) a«rf 

u(^-V/ 2 -a(|ur™~ 1 )-V^)da;dt + J u o (x)ip(0,x)dx = (5) 

/or eacft test function tp G C 1 (Qt) swc/i iftai </? /las compact support in the space variable x, and vanishes 
near t = T . 

Theorem 2.2 (Self-similar solutions) For each a G (0, 2], m > 1 and i? > 0, i/ie function 

® am (y)=(k(R 2 -\y\ 2 )l)^ with fc=f 3; ^ )( -, , jt A (6) 

' v +/ V^-i) + a/ V 2Qr ( 1 + f) r (^)/ 

is a solution of (4). Consequently, the function 

u(t, x) = £-<*(".-!)+<* ( A: (i? 2 - |x| 2 t""<"- 2 i)+°) 2 ) (7) 



is a weafc solution of (1), satisfies the equation in the pointwise sense for \x\ ^ Rt d <.™- 1 ~>+ a , and i 



is 



min | 2(m-i) ' 1^ -Holder continuous at the interface \x\ = i?i d ( m - 1 '+° . 

When a = 2, wc recover the classical Kompanccts-Zcl'dovich-Barenblatt-Pattle formulas, see, e.g., [13]. 
Note also that given M > there exists a unique R > such that J" ^>a,m(y) dy = M. 

Remark 2.3 As mentioned above, self-similar solutions of (1) have been proved to exist in [3] by studying 
the following obstacle problem for the fractional Laplacian 

P>$, V=(-A)$P>0, either P = $ or V = 0, 

itrei/i a G (0, 2) and $(2/) = C — a|j/| 2 . The novelty of our approach is that we exhibit the explicit solution of 
this obstacle problem: P(y) — I a ($0.2) (y), where I a is the Riesz potential, and $ Qj2 is defined in (6) with 
m = 2 and a suitable R > 0. 

Those explicit self-similar solutions express one of the most important features of the porous medium 
equation: the property of finite propagation speed. In the case of the classical porous medium equation 
(q = 2), this property is established using comparison with suitably large self-similar solutions, cf. [13]. 
For the generalized porous medium equation (1) with m = 2, special supersolutions have been used for 
comparison, cf. [4]. 

The proof of Theorem 2.2 is based on an application of the following fundamental technical fact. 
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Lemma 2.4 For all /3 G (0,2) and 7 > 0, we have 

„ \c jAd , 2 Fj d ^,-l±l-, d -- A yA br|„|<l, 

I C 7 ,/3,rf|yl x 2^1 I — jj— .—3— Ir-ja ) for [y| > 1, 

r(x+i)r(£=£) ~ r(~* \ i)r( d ~ ? ) 

with Cy^d = 2 p7|WgE^"jl ar ^ ^r^.d = 2 -/3 r ^ r ^d+7 , a » w/iere %Fi i^e hypergeometric function. 

The verification of (8) consists in passing to Fourier transforms and calculating certain integrals (the so 
called (Sonine-)Weber-Schafheitlin integrals) involving Bessel and hypergeometric functions 2 F\, cf. [11]. 

Proof of Theorem 2.2. Let <f> a (y) = (l - |y| 2 )|- Observe that ^ € L 1 (R d ). We next show that h- a (<t>a) G 
Wioc( Rd )- B y Lemma 2 - 4 with 7 = Qg(0,2), J 3 = 2-a, and 2 .Fi(a, -1; c;z) = l- §z, we get 



if |tf|<l, 




h-a((/>a)(v) = 



which is a locally intcgrable function. Recalling that V" -1 = Via-tm we then deduce by the chain rule that 
for y G B\, V a ^ 1 ((j> a )(y) = — (dKa^) 1 y where K a> d is defined in Corollary 2.5 below. For \y\ > 1, one 
uses J| 2-^1(0, 6; c; z) = 2 F\(a + 1, 6 + 1; c + 1; z), hence V Q_1 (0 Q ) is locally integrable. To conclude, we 
remark that 

0S^ T (tf)V a -Va(») = -^{V) (Jfa.dd) -1 V, for all ?/ G M d , 
because <^ Q (y) = for |y| > 1. 

Hence, scaling the variables, we immediately obtain that the function § a ^ m defined in (6) satisfies 
ya-i^m-i) = for | y | < R NoWi for all y g R d the identity 

follows because $ a , m (y) = for all \y\ > R, so (4) holds with A = d ( m _\) +a and k = dX K a ^d defined in (6). 

It is straightforward to verify using (4) that u given by formula (7) is a weak solution of (1) in each strip 
(to, T) x M. d , < to < T < 00. Moreover, the Holder continuity is easy to check. □ 

The following known result (with an important probabilistic interpretation) proved by Getoor [6, Th. 
5.2], see also [9, App.] for a related calculation, is an immediate consequence of Lemma 2.4. For its proof, it 
suffices to use the relation (— A)t = AI 2 - a - 

Corollary 2.5 For each a G (0,2], the identity K a: d(—/S)^ Ml — \y\ 2 ) _? ) = —1 in B\ holds with the 



r (-) 

explicit constant K a>d = 2 ~ r (i 1 °)r( d+ °) ' 



3. The Cauchy problem and asymptotics 



We now briefly discuss the questions of the existence of weak solutions and their uniqueness in the case 
m = 2. In [1], viscosity solutions have been considered which permitted the authors to prove regularity and 
uniqueness of solutions to (1) in one space dimension. In higher dimensions, a construction of mild solutions 
is achieved in [3] through a parabolic regularization of (1) with a suitable cutoff of the singular kernel of 
li-a, and then a passage to the limit. The uniqueness of weak solutions and, a fortiori, the validity of 
the full comparison principle seem to be difficult questions, cf. a discussion in [3]. Another construction of 
weak solutions to (1) with m > 1 as limits of mild solutions of parabolically perturbed equation (1) will be 
published in [2]. 
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Formal computation 3.1 (Decay of solutions for the Cauchy problem) Suppose that for < «o G 
L 1 (R d )nL°° (R d ) . Then L p (M. d ) norms (1 < p < oo) of any sufficiently regular and global in time nonnegative 
weak solution u of (l)-(2) such that J u(t,x)dx = J Uo(x)dx decay algebraically 

n + d(m-l)/p d 

\\u(t)\\ p < C(d, a, r7T,,p)||wo|| 1 ° +!i<m " 1) t~ J i"- 1 )+° (1 ' 1/p) forallt>0. (9) 



Our computation below can be applied, for example, to weak solutions constructed by Caffarclli and 
Vazquez [3], who studied problem (l)-(2) with m = 2 and with initial conditions satisfying < uq(x) < 
Ce~ c \ x \ for some constants c, C and all x e R d . 

Formal proof of (9). In order to prove the announced LP estimates of solutions (similar to those for degen- 
erated partial differential equations like the porous medium equation in e.g. [14], [5, Ch. 2]), we recall the 
Stroock-Varopoulos inequality for q > 1 

v\ q - 2 w(-A)^wdx > 4( - g ~ 2 V> J V^|w|5 2 da; (10) 

valid for each w £ L q (R d ) such that (— A)%w <G L q (M. d ). Note that the constant in (10) is the same as for 
the usual Laplacian operator —A (i.e. a = 2). The proof is given, e.g., in [10, Prop. 1.6] and [10, Th. 2.1 
combined with (1.7)]. 

We will also need the following Nash inequality 

\Hl il+ ^ <C N \\V%v\\j\\v\\f (11) 

valid for all functions v with v e L 1 (M d ), E L 2 (R d ) with a constant C N = C(d,a). The proof of (11) 

for d = 1 can be found in, e.g., [8, Lemma 2.2], and this extends easily to the general case d > 1. 
Moreover, we will need the Gagliardo-Nirenberg type inequality: for p > 1 and p > m — 1, 

IW<Civ||V*|«|5||a|H|S (12) 

with a = ~f d ( r ~p +a and b = ^"^^y 1 ^ and r = p + m — 1. This inequality is a consequence of the Nash 

r II a r ||2 ra 

inequality (11) written for v = , i.e. \\u\\ r d < Cn ||^ t I u I^|| 2 IMIr , and two Holder inequalities for 
the L* norms: ||t*|| p < HI?INi~ 7 and \\u\\r < |M|£|M|J- 5 , which hold with 7 = ^ and S = r=f 2. 
Combining the above three inequalities, we get (12). 

The computation, which lead to (9), consists in getting a differential inequality of the form -4r J \u\ p dx < 
—K\\u\\p\\u\\^ b for some positive constant K and where a and 6 appear in (12). 

Multiply (1) by u p_1 with p > 1, integrate by parts, and use the relation V • V Q_1 = — (— A)^ to get 



1 d 
p~dt 



/•^--(p-d/^v-v-V-')* 

= - V — - [ u p (—A)^u m ~ 1 dx 
P J 

4(p-l)(m-l) „a / P+ m -l N 2 
< - V~ 



(p + m — l) 2 
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after applying the Stroock-Varopoulos inequality (10) with w = u m 1 and q = — ^ + 1. To estimate the 
right hand side of the above inequality, we use (12). We finally have: 



A J \u\ p dx<-K\\u\\;\\u\\r 



with K = 



_ 1 4p(p-l)(m-l) 



The above inequality leads to the differential inequality -§^f{t) < —KM b f(t)p 
for the function f(t) = \\u(t)\\P, M = ||wo||i, and a/p > 1, which immediately gives the algebraic decay of 



C N (p+m-l) 2 



the LP norms for p>m-l: f(t) < [K ( - — 1) M~ b t 



Finally, we obtain the desired estimate with 



C(d, a, m,p) 



K 



4(to - l)(d(m - 1) + a) 



(p + m — l) 2 



Of course, this is sufficient to get the conclusion of Theorem 3.1 for each 1 < p < oo since ~~ 



and interpolating between L 1 (K d ) and L p (M. d ) with p sufficiently large. 



d(m-l)+a' 

□ 



Remark 3.2 In our recent work [2], we present a more subtle iterative argument, which allows us to show 
that also 

IKOIloo < (7|Klli (m " 1,+Q i~ d(m - 1)+£ * for all t > 0. 
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